Self-Reduction Rate of a Microtubule 
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We formulate and study a quantum field theory of a microtubule, a basic element of living cells. 
Following the quantum theory of consciousness by Hameroff and Penrose, we let the system to reduce 
to one of the classical states without measurement if certain conditions are satisfied(self-reductions), 
and calculate the self-reduction time tat (the mean interval between two successive self-reductions) 
of a cluster consisting of more than A'' neighboring tubulins (basic units composing a microtubule). 
tm is interpreted there as an instance of the stream of consciousness. We analyze the dependence 
of rjv upon A'' and the initial conditions, etc. For relatively large electron hopping amplitude, rjv 
obeys a power law rjv ~ A''*', which can be explained by the percolation theory. For sufficiently 
small values of the electron hopping amplitude, rjv obeys an exponential law, rjv ^ exp(c'A'^). By 
using this law, we estimate the condition for rjv to take realistic values rjv ^ sec as A'' ^ 1000. 



I. INTRODUCTION 

Most approaches to understanding various character- 
istic functions of the human brain in a framework of con- 
ventional physics are classical. For example, the Hopfield 
model of neural network [1], a standard model of asso- 
ciative memory, adopts an Ising spin variable (±1) to 
describe the state of each neuron. In the field of psy- 
chophysics, phenomenological relations between psycho- 
logical quantities (like the sense of pain) and physical 
quantities (like voltage) in classical physics are of main 
concern. Although most researchers may agree that the 
human brain should be described by quantum theory of 
electrons and atoms at the microscopic level, there seems 
to be a prejudice that the functions of human brain in 
our daily life are to be well understood, if possible, by 
notion of classical physics without recourse to quantum 
physics. 

The problem how to understand mind and conscious- 
ness remains as a central and still open question in the 
physics of the brain. We don't know even whether the an- 
swer to this problem exists. In such a situation, Hameroff 
and Penrose [2] proposed a quantum theory of conscious- 
ness according to a new quantum mechanics of micro- 
tubule. Their idea is interesting because it suggests the 
fundamental importance of quantum effects like quan- 
tum measurements and collapses of wave function upon 
the functions of human brain. 

Microtubules are basic building blocks of living cells 
including neurons and axons, conveying mitochondria, 
etc. Each microtubule is a hollow cylinder with a di- 
ameter ~ 25nm, the surface of which is a 2D array of 
units called tubulins (See Fig. la). The array consists of 
13 columns, and its longitudinal size L is not definite 
(typically L = 100 ^ 1000). Neighboring tubulins in two 
adjacent columns is displaced with a fixed pitch, so that 
tubulins form a triangular lattice. A tubulin is a dielec- 
tric dimer consisting of an a tubulin monomer and a (3 
tubulin monomer. 



As a quantum model of a tubulin it is often treated 
[2-5] as each tubulin contains a single active(mobile) elec- 
tron and takes two independent states, \a) and ac- 
cording to the "upper" or "lower" position of the elec- 
tron (See Fig. lb). Tuszynski et al. [3] set up and studied 
some quantum models of a microtubule as a conventional 
quantum system, i.e., without self- reductions. Also, Ras- 
mussen et al. [4] stuidied signal propagations along a mi- 
crotubule by using a network model of cellular automata. 

Penrose [5] argued that any quantum system is to be 
affected through couplings to quantum gravity in such a 
way that its wave function should make self-reductions 
(called orchestrated objective reductions) according to 
certain rules that reflect the uncertainty principle in 
quantum gravity. In the theory of Hameroff and Penrose 
[2], a well developed coherent (i.e., superposed) state of 
each tubulin, Ca\a) -\- Cp\P) with \Ca\ ~ \Cp\ has a too 
large uncertainty of the location of electron, and may be 
unstable if some other tubulins are also in such states at 
the same time, from the viewpoint of uncertainty princi- 
ple of quantum gravity. So each of these tubulins should 
make a self-reduction to be put back to one of the eigen- 
states of electron position, \a) or |/3). They argue that a 
sequence of such reductions works as a clock and forms a 
stream of consciousness. The mean time between two 
successive reductions, which we call the self-reduction 
time r, is interpreted as each moment of consciousness. 

Then quantitative estimation of r would become a 
main concern. However, Tegmark [6] posed a ques- 
tion that a microtubule suffers from decoherence effects 
through quantum interactions with its environment, for 
instance ions and water molecules in the surrounding cy- 
toplasm, and the decoherence time Tdcc (the time scale 
to lose off-diagonal elements of the reduced density ma- 
trix, hence to destroy long-range quantum superposi- 
tions) may be much smaller than the self-reduction time 
T. If this is true, the scenario by Hameroff and Penrose 
loses its reality. He estimated the decoherence time Tdec 
of an entire microtubule at finite temperatures as rjcc ~ 
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FIG. 1. (a) Development of a microtubule where each sec- 
tor represents a tubulin and (b) Orientation of the six neigh- 
boring tubulins (from Ref.[3]; lengths in unit of nm). Each 
filled circle represents a possible location of electron with cor- 
responding quantum state \a) (upper) or |/3) (lower).® 

10~^^ sec and claimed that microtubules can be treated 
classically in the interesting time scales 10~^ ^ 10~^ sec 
for a stream of consciousness, and quantum effects play 
no significant role there. Then Hagan et al. [7] responded 
to Tegmark by arguing in details that his estimate is 
not definitive; Tdoc may be much larger and the quantum 
theory may be still relevant. 

In this paper we leave this problem of estimating the 
decoherence time as a still controversial and open prob- 
lem, and concentrate on the quantitative estimation of 
the self-reduction time r itself of an isolated microtubule 
at zero temperature without external disturbances, thus 
ignoring decoherence effect by the environment. An en- 
semble of microtubules and/or the effects of temperature 
and environment [6,7] should be considered in the next 
step. Explicitly, we set up a quantum theory based on 
the classical model of Rasmussen et al. [4] and calculate 
T by proposing some rules of self-reductions based on the 
Hameroff-Penrose theory. Such a study should certainly 
present some important informations to scrutinize the 
relevance of Hameroff-Penrose theory. Since the typical 
time scale set by the Coulomb energy at nano scales is 
^ 10~^^ sec, the central question is whether r can take 
values of the order of 10^^ ^ 10^^ sec which seems to 
be reasonable as the moment of consciousness. We shall 
see that r exhibits systematic behaviors which may be 
typical for models of self-reductions. 

One may conceive reasons of self-reductions other than 
Penrose's proposal based on quantum gravity [5]. Ghi- 
rard et al. [8] propose a modification of Schrodinger 
equation to a nonlinear stochastic equation to incorpo- 
rate dynamical reductions. As yet another possibility, 
(quasi)reductions may occur due to internal quantum 
measurements of microtubule by surrounding environ- 
ment in the brain, or even the microtubule itself may 
"measure" its part [9]. Our rules of self-reductions are 
phenomenological and may mimic these cases. 

The rest of the paper is organized as follows: In Sect. 2 
we set up a quantum field theory of a microtubule. It 



may be interpreted as a frustrated quantum spin model 
in two dimensions. We propose to study its dynamics by 
solving an approximate Schrodinger equation derived by 
the variational method. In Sect. 3 we calculate the self- 
reduction time Tjv for various conditions. We find that 
tn exhibits different behaviors depending whether the 
electron hopping amplitude is large or small. In Sect. 4 
we present conclusion. 



II. MODEL 



A. Hamiltonian 



The Hamiltonian H of a microtubule is given by 

(« j) 7,7' 



7 — Q,^ 



(2.1) 



where a = (3, (3 = a. bij, b]^ (7 = a, (3) are the fermionic 
annihilation and creation operators of the electron in the 
7 state at the i-th tubulin, satisfying 



1^17, ^j:'7']+ — 0, [6i7, — (5ii'5-y^'. (2-2) 

e is a dielectric constant of the microtubule. The first 
term of H represents the kinetic energy of the electron 
changing its state a [3 within the i-th tubulin with the 
hopping amplitude k. The second term represents the 
Coulomb energy V^^ between two electrons at (17) and 
(j7') separated by the distance i?^/ , the value of which 
is given in Ref. [4]. j) implies the sum over all the 
neighboring pairs i,j given in Fig. lb (Each tubulin has 
six neighbors). The number of electrons in each tubulin 
is a constant of motion, 



N, = Y,b\^b,^, [i?,iV.]-0. 



(2.3) 



In actual calculations, we start and stay in the subspace 

N^ = 1. 

This Hamiltonian may be expressed in terms of 
the s = 1/2 SU(2) quantum spin operator Si = 
{h/2){bl^,blp)a {bia,bipf {a are Pauh matrices) as 



H = 



/ y JjjSjzSjz ^ ^ BSj 



const. 



(2.4) 



with 



J,. = - 



V, 



v. 
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The values of J^- are 

■/North = ■/south = 0.0833?l~^ Vq , 

JsE = Jnw = 0.0091ft- Vo, 
Jne = Jsw = -O.G280fi-Vo, 



Vo 



1 



Inm 



• 2.31 X IQ-^^ Joule, 



(2.6) 



where Vo is the Coulomb energy of a pair of electrons 

in the microtubule separated by 1 nm. So the system 
may be viewed as a frustrated spin model in an external 
magnetic field B. At fc = 0, the system involves only 
^2 , and its groimd state is foimd to be the "stripe state" 
in which Sz are aligned to ±fi,/2 along each column with 
the alternative signs, except for a pair of two NN columns 
(1st and 13th, say) with the same signs. The degeneracy 
is 13 X 2 [10]. As k increases, the Sx component develops, 
and at — > cx), all Si align to (/i/2)(l,0,0). 



B. Time evolution in variational method 

In the time interval between two successive self- 
reductions, the state vector of ^ microtubule at 
time t evolves according to the Schrodinger equation, 



(2.1) 



If the microtubule contains V tubulins, then \tp{t)) is 2^ 
dimensional, and for F ^ 13 x 100, the precise evaluation 
of \ip{t)) is beyond the ability of our computers. Thus, we 
evaluate |f/'(0) approximately by the variational method. 
Explicitly, we choose the variational state \il>v{t)) in the 
factorized form, 

\Mt)) = n [ '^-(*) \ic^)+Ci0{t) m ] , 

i 

\Cia{t)\^ + \Cii3{t)\^ = I, (2.2) 



where 



(2.3) 



is the state of the i-th tubulin in which the electron is in 
the 7 state. By minimizing the action, 

S = j'j dt (V„(OI (^hJ^ - \Mt)), (2.4) 

we get the following equations of motion for 2V complex 

iicic 



coefficients Ci^: 



ih- 



dt 



Eq.(2.5) respects the unitarity, 

\Cic,{t)\^ + \Cip{t)\'' = l, 



(2.6) 



and keeps the average Ey = {tjjy{t)\H\tpy{t)) a constant 
of motion. 

We note that, if AF-] = F^]" - V^"/ is negligible, 

^^7^jE(^J" + ^if) = ^~ (2-7) 



becomes independent of neighboring Cjj's, so Eq.(2.5) 
are decoupled to give the mean- field solution (MFS), 



a. 



Ci+e-'^+' + (-)Ci_e-'"-*, 



Ci± = - [Ci40) ± Ci0{O)] , n±=JTk. (2.8) 

The actual value of |AT/J|/(V;]" + yj^) averaged over 
j is 0.19, so the couplings between neighboring tubulins 
are intermediate. 

We solve the differential equations (2.5) numerically 
by the Runge-Kutta method with discrete time step At. 
As the initial condition Ci^{0) at t = 0, we consider the 
following four cases: 

(i) US(Uniform Start); C,a(0) = 1,Q^ = for all i, 

(ii) SS(Stripe Start); Cja(O) = 1,0 for alternative 
columns as explained. 

(iii) RIS(Random Ising Start); Ciy(O) — 1, 0, randomly, 

(iv) RCS (Random complex Start); Ci-y(O) G C randomly. 



C. Rule of self- reductions 

For the rule of self-reductions, we follow the idea of 

Hamcroff and Penrose [2]. We first choose a parameter 
Pq and prepare the "reset zone" of Cia(t), 



Po < \ac.{t)\^ < Pi{= I - Po) 



(2.9) 



j (NN to i) 



When Cia(t) is in this zone, we judge that the state \ia) 
is sufficiently coherent and may need a self reduction. In 
the simulation we watch whether |CM(t)P falls into this 
zone at every time step. The tubulins that arc in the reset 
zone form a cluster, a set of tubulins that arc connected 
each other as neighboring pairs. If the size of cluster at 
t = is same or more than the prefixed number A'', we 
let all tubulins within this cluster make self-reductions at 
the next time step. As the rule to determine Cj-y(ti{+At), 
we consider the following two cases: 

LR : Local reduction; Each tubulin reduces to \ia) or \ifi) 
according to the probabilities \Ci^{tR)^. 

GR : Global reduction; All the tubulins in the cluster reduce 
to the common state, I7), with the probabilities 
{\Ci'y{tB)\^), the averages over the cluster. 
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tCxlO'^At) 8 4 t(xlO=At) 
(a) (b) 
FIG. 2. Typical behavior of \C\a\'^{t) in a run for RCS/LR, 
L = 100. (a) k = O.ll^cA^ = 100 (b) k = 0.01Vo,Af = 30. 

In Fig. 2 we show typical time dependence of \Cia{t)\^ 
for two values of k. For larger k the smooth motion be- 
tween reductions is rather regular, reaching the endpoints 
and 1 at every period. For smaller fc, the motion be- 
tween reductions is often pulled back before reaching the 
endpoints. This difference is explained as the effect of 
potential term in H (See Sect. Ilia and Illb). 

In Fig. 3 we show typical snapshots of \Cia{t)\'^ for a 
part of a microtubule before and after a reduction. 

To calculate the self-reduction time tn for the mini- 
mum cluster size iV, we make a sufficiently long run with 
the period ^totai for each N. is then defined as 



TN 



lim 



^total 



(2.10) 



where is the total number of all the self-reductions 
within the run. So tjv is the average time interval be- 
tween two successive reductions, tn may be expressed 
by using the probability pN that clusters of size equal to 
or larger than N appear as 



\m—l ^ 



7n—l 



1 

Pn' 



(2.11) 



In the simulations we consider microtubules of the 
length L up to L = 2000 (16 /im) with the open boundary 
condition. Wc choose 




At = 0.01 X to, 

_ h he y. Inm 
*o = t7 ^ 2 



4.571 X 10"'^ X e sec, (2.12) 



where to is the time scale set by the Coulomb energy Vb 
at 1 nm. We choose sufficiently large itotai up to 



itotai = 2 X lO^steps X At, 



(2.13) 



to obtain stable values of rjy. The above At is checked 
to be sufficiently small for this ttotai- To determine the 
clusters at each t, we adopt the algorithm of Hoshen and 
Kopelman [11] developed in percolation theory [12]. In 
Table 1 we collect the parameters and their values used 
in our simulations. 



parameter 


symbol 


value 


hopping 


k 


0.002 - 2(xVb) 


amplitude 






length 


L 


50 - 2000 




miiform (US) 




initial 


stripe (SS) 


All 


condition 


random Ising (RIS) 






random complex (RCS) 




cluster size 


N 


5 - 26000 


reset zone 


Po 


0.01 - 0.49 


reduction 


global reduction (GR) 




condition 


local reduction (LR) 


both 


elapsed time ttotai 


10^ - 2 x lO'^(xto) 



(b) t=506At (c) t=584At (d) t=585At 



Table 1. Parameters used in simulation. 

We note that the Schrodinger equation (2.1) is invari- 
ant under t ^ At and H — *■ X^^H. This is refiected also 
in the equation of motion (2.5). Thus, r is proportional 
to the dielectric constant e when k is measured in unit of 
Vq. Experimental measurements of e are still not defini- 
tive but predict values like e = 8.41 [13], 1 < e < 100 
[14], etc. Also, if the charge of the mobile object which 
determines the tubulin state, \a) or \f3), is modified from 
— e of a single electron to —ne, its effect is to be reflected 
by replacing e ^ e/n^ in the formulae below. 



III. SELF-REDUCTION TIME 

As we shall see, exhibits different N dependences 
for large fc's and small fc's. We study the case of large k 
first, and then the case of small k. 



FIG. 3. Snapshots of time evolution of |Cia|^ in a part of 
a microtubule. Darkness of each tubulin represents jCiaj'^ 
between O(white) and 1 (black) in gray scales. (a) RIS 
at t = 0, (b) coherent evolution, (c) sufficient number of 
tubulins are just in the reset zone (a cluster is formed), 
(d) most of tubulins are reset by the GR to |/3) (white). 
k = O.lVb.A'' = 200, L = 100. 



A. Large hopping amplitude k 

Let us first consider the case of large k. In Fig. 4 we 
plot ttv vs. N for k = 0.1 Vb for several values of Po- For 
all the four initial conditions, tat up to can be 

fitted well by the power law. 
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Tsync (Po=0.4) : 
Xsync(Po=0.1)- 




* Po=0.4 

A Po=0.3 

X Po=0.2 

O Po=0.1 



FIG. 4. TN VS. AT for fc = O.lVo and L = 500 with RIS/LR. 
Errors are estimated by 20 independent runs. The dashed 

straight lines show the power law of (3.1). The horizontal 
lines show the saturation value Tsync defined by Eq.(3.7). 



rjv — aN'' X to, 
a ~ 1.02 X 10-%, b ~ 1.08 for Pq = 0.3. 



(3.1) 



Eq.(3.1) may be explained by adopting the independent 
oscillator model where \Cia{t)\^ is given by (2.8). Then 
the probability p that each |C8a(t)|^ is in the reset zone 
is given by 



p = l--cos-^(l-2Po)- 

TT 



(3.2) 



For Pq > 0.146, p is below the critical probability = 0.5 
of the triangular lattice above which percolation takes 
place [12]. The scaling argument [12] predicts pn for 
p<Pc as 



p^oc^^s /(z)oc|^_,^^ {p^-p)N^>l, ^^-^^ 

s—N ^ 

where f{z),z = [pc ~ p)s'^, is the scaling function and 
the exponents are given by [12] 



187 36 

— , a=-, coc(p,-p)/ . 



(3.4) 
(3.5) 



By assuming 

{pc-p)N^<l, 

tn oc 1/pn oc N^~^ giving us the form of (3.1) with 

6 = 5- 1 = 1.055. (3.6) 

For only one large cluster is possible and all 

the tubulins in it exhibit a synchronized (collective) mo- 
tion as 



\Cin 



k k 
Psync(t) = cos^[-(t - ffl)] or sin^[^(* - tn)], 



TN ^ Tsync = — X cos ^(1 - 2Po), 



(3.7) 



as the MFS (2.8) predicts. In fact, Fig.4 shows that tn 

saturates to r^ync for N>V/2. Thus k = 0.1 Vq is "large", 
where the coupling to neighbors via AV^J is negligible. 
Here the system is in the single-body regime. 



B. Small hopping amplitude k 



Let us next consider the case of small fc's. For small 
fc's, Al^J becomes relevant. In Fig. 5 we plot tn vs. A'' at 
k = O.OIVq. Here each |C,c«(t)|^ basically sweeps between 
and 1 with a basic frequency ^ k/h, but is sometimes 
"pulled back" due to AP"J (Sec Fig.2b). The system is 
in the many-body regime where the potential term is rel- 
evant. As A'' increases, tn in Fig. 5 exhibits a crossover 
from the power law to the exponential law. We note that, 
although the tubulins arc not independent each other 
here as explained, Eq.(3.3) may describe this crossover 
phenomenologically (with an effective Pc — p, etc.). The 
exponential behavior of Fig.5 is fitted as 



Tjv — a' exp{c'N), 
a' ~ (0.40 ~ 1.30) X 10"" x e sec, 
c' ~ 0.025 ~ 0.029 for Pq = 0.3, L = 



50 ^ 200. (3.8) 



To study the condition for tjv to take realistic values, 
extremely long runs are required. So let us assume (3.8) 
to hold for larger iV's as an asymptotic expression, and 
extrapolate it to estimate the minimum size N of the 
cluster. For /c = 0.01^0, Po = 0.3 and n' = 0.85 x lO"!"* x 
e sec, c' = 0.027 we obtain the minimum N for typical 
values of the dielectric constant e in Table 2 below. 





T > 10"^ sec 


T > 10"^ sec 


e = 1 


N > 1030 


N > 1115 


e = 10 


N > 945 


N > 1030 


e= 100 


N > 860 


N > 945 



Table 2. Minimum size of the cluster for r > 10 ^sec, 
10~^sec and various e at fc = O.OlVb and Pq = 0.3 calcu- 
lated by Eq.(3.8). 



L=50, 100,200 




100 200 N 

FIG. 5. TN vs. AT for fe = 0.01 Vo with RIS/LR. tn exhibits 

a crossover from the power law (oc A^*" , b' ~ 1.56 for Pq ~ 0.3; 
dashed curves) to the exponential law (3.8) (dashed straight 
lines). 
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0.01 0.1 1 (xVO) 



FIG. 6. TN vs. fc for L = 200 and Pq = 0.3. Upper data in 
black are for N = V and lower data in gray are for A'^ = O.iy. 
Dashed line is rgync of (3.7). 




0.1 0.2 0.3 0.4 0.5 

Po 



FIG. 8. Two regimes in the Pq — k plane separated by 
k — ka{Po)- The smooth curve separating two regimes is 
an interpolation of fc = ka{Po) determined from Fig. 7 (six 
small circles). 



C. fc-dependence and the estimation of k 

Let us study the fc-dependence of t^v in detail. In 
Fig. 6 we plot Tjv vs. fc for N = O.IV and N = V with 
L = 200, Po = 0.3 and LR. Let us first consider N ^V. 
Simulations with RCS show no reductions throughout 
runs due to the random initial configurations. For the 
other three initial conditions, the data (in black) at large 



fc's lie on the line rjv 



(3.7) of the synchronized 



behavior of single cluster explained for Fig. 4. As fc de- 
creases, they deviate from Tgync suddenly at the point 
marked as a with k = ka — 0.07Vb As fc ^ 0, tjv blows 
up. This is because these SS, US, and RIS become ex- 
act eigenstates of at fc = 0; each initial eigenstate has 
time-independent \Cia{t)\'^{= 0,1), hence no reductions 
are possible. Let us next consider N = O.IV. The SS, 
US, and RIS data start from Tgync at large fc and devi- 
ate from it downward at the point b and join the RCS 
data. They start to blow up at c. The decrease of tn at 
b reflects generations of plural clusters. 

In Fig. 7 we plot r^r vs. fc for various values of Po with 
RIS/LR, N = V/2, L = 200. The data show systematic 




]^ l(xvo) 

FIG. 7. TN vs. k for various values of Po, i = 200, TV = V/2 
and RIS/LR. Dashed lines are Tsync of (3.7). 



dependence on Pq. As Pq increases, the reset zone be- 
comes narrow and tj\j at fixed fc increases. The point of 
blowing up, fc = fca, also increases as Pq. 

In Fig. 8 we plot an interpolating curve of fca in Fig. 7 
as a function of Po. It separates two regimes: (i) single- 
body regime for k > ka where the kinetic term in H 
dominates over the potential term and the dynamics of 
each tubulin is well described by the synchronized motion 
of Eq.(3.7) and (ii) many-body regime for fc < fc^ where 
the potential term dominates over the kinetic term. As 
Figs. 6, 7 show, the change between these two regimes is 
very sharp. 

Let us estimate the value of fc. For this purpose, we as- 
sume a hydrogen-like wave function for \a) and |/3) states; 



(r|7) cx exp( 



e 



(3.9) 



Then, by a straightforward calculation, one obtains 
, (a|^|/3)-(a|/3)(a|^|a) 



6 (1 - {am 



(3.10) 



In Fig. 9 we present fc vs. £ calculated from Eq.(3.10) 
by putting — r^|=4nm as shown in Fig. lb. It seems 
that the reasonable value of fc may be within the range 
O.OlVb < fc < O.lVb. 




I L.D z 2.5 (nm) 

FIG. 9. k vs. £ given by Eq.(3.10). Numerically, it gives 
k = O.lVo for £ = 2.0nm and k = O.OlVo for £ = 0.61nm. 
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IV. CONCLUSION AND DISCUSSION 

If we write the time period of a physiologically reason- 
able value of each moment in the stream of conscious- 
ness as Tconj a necessary condition that the scenario of 
Hameroff and Penrose [2] becomes successful reads 

"^dec ~ Tat ~ ''con- 

(4.1) 

For 

'^con it is reasonable to estimate as Tcon 
IQ-^ sec. Then the exponential law (3.8) for tat at small 
k shows that, in order to satisfy the second inequality 
of (4.1), one needs N > 1000 as Table 2 shows. This 
cluster size is realistic enough, because it is the same 
order or smaller than the typical size of microtubule, 
V = (100 ~ 1000) X 13, explained in Sect.l. Thus, our 
present analysis shows that the second inequality holds. 
Concerning to the first inequality of (4.1), estimation of 
Tdec is crucial, but, as stated in Sect.l, it is still contro- 
versial [6,7]. We need reliable estimation of Tdec- Here 
we point out that a possible quantum ordered state of 
a microtubule or a set of microtubules with spontaneous 
symmetry breaking and an off-diagonal long-range order 
may give rise to an extra stability of coherence against 
decohcrent fluctuations due to the environment, leading 
to much larger values of Tdec- This possibility is pointed 
out by Penrose [5] and others in a general point of view. 
Because there is now an explicit quantum model (2.1), 
one may explore such possibility in a quantitative man- 
ner using conventional techniques of quantum field theory 
and statistical mechanics. This is certainly an interesting 
future problem. 

In conclusion, we set up a quantum field theory of 
microtubule and estimated the self-reduction time rjv- 
Tjv exhibits systematic behaviors; the power law in the 
single-body regime for k > ka, or the exponential law 
in the many-body regime for k < ka- tn can take ar- 
bitrary large values if one imposes unrealistic conditions 
like k 0, A'^ oo, or Pq — > 0.5. Further experimen- 
tal studies of microtubules, e.g., determination of k, may 
lead us to better understanding of their behaviors; We 
may then be ready to answer some interesting questions 
like whether they are sitting on "the edge of chaos" dis- 
cussed in the study of complex systems [15]. 

Concerning to the approximation(variational method), 
we point out a possibility that the approximate solution 
plus self-reductions as we imposed may mimic the exact 
solutions of Schrodinger equation. Then we do not need 
extra coupling to quantum gravity claimed by Penrose. 
It is interesting to investigate whether the conventional 
quantum theory itself exhibits (quasi)reductions by let- 
ting its subsystems make sclf-(quasi)reductions [9]. 
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